function [mid_T,uend,u,edge_T]=PDE(alpha,h,u0,T0,time)
    d=0.15*1e-3;dx=d/2;%芯片温度计算中心及两边界
    dt=diff(time);dt=dt(1);
    r=alpha^2*dt/(dx^2);
    n=length(time);
    u=zeros(3,n);
    u(:,1)=u0;
    for j=1:n-1
        u(2,j+1)=r*u(1,j)+(1-2*r)*u(2,j)+r*u(3,j);%热传导方程计算中心温度
        u(1,j+1)=(u(2,j+1)+h*dx*T0(j+1))/(1+h*dx);%边界条件计算边界温度
        u(end,j+1)=u(1,j+1);%对称性      
    end
 
    mid_T=u(2,:);%中心温度
    mid_T=mid_T';
    uend=u(:,end);%下段初值
    edge_T=u(1,:);%边缘温度
end
